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ABSTRACT 

We present a population synthesis study of NGC 4435, an early- type Virgo galaxy interacting with 
NGC 4438. We combine new spectroscopic observations obtained with the Spitzer Space Telescope 
IRS instrument with IRAC archival data and broad band data from the literature. The IRS spectrum 
shows prominent PAH features, low ionization emission lines and H2 rotational lines arising from the 
dusty circumnuclear disk characterizing this galaxy. The central SED, from X-ray to radio, is well 
fitted by a model of an exponential burst superimposed on an old simple stellar population. From the 
lack of high excitation nebular lines, the [Nein]15.5/[Neii]12.8 ratio, the temperature of molecular 
hydrogen, and the fit to the full X-ray to radio SED we argue that the present activity of the galaxy 
is driven by star formation alone. The AGN contribution to the ionizing flux is constrained to be 
less than 2%. The age of the burst is found to be around 190 Myr and it is fully consistent with the 
notion that the star formation process has been triggered by the interaction with NGC 4438. The 
mass involved in the rejuvenation episode turns out to be less than 1.5% of the stellar galaxy mass 
sampled in a 5" central aperture. This is enough to render NGC 4435 closely similar to a typical 
interacting early- type galaxy with inverted Can[H+K] lines that will later turn into a typical cluster 
E+A galaxy and enforces the notion that these objects are the result of a recent rejuvenation episode 
rather than a genuine delayed formation. 

Subject headings: galaxies: elliptical — galaxies: evolution — galaxies: stellar content — galaxies: 
individual (NGC 4435) 



1. INTRODUCTION 

NGC 4435 (VCC1030), a n early-type SB0i(7) galaxy 
l)Sandage & Tammanlll98'Tj) member of the Virgo Clus- 
ter, is of great interest because it is a nearby exam- 
ple of galaxy interactions in the cluster environment. 
NGC 4435 is experiencing an off-center encounter with 
NGC 4438, a bulge dominated late-type galaxy lying at 
a projected separation of 4.3' (i.e. 20 kpc, adopt ing a 
Virgo distance of 16.1 Mpc from Kel son et al.l2000]) . Ev- 
idence of the strong on-going interaction comes from the 
well developed stellar tidal tail («30 kpc) of NGC 4438 
but also from the structure of the inter s tellar medium 
(ISM ) of both galaxies. Ha+[Nn] images l|Kennev et al.l 
1995) reveal several ionized gas filaments starting from 
the disk of NGC 4438 and extending out of the disk up 
to 10 kpc. The study of the distribution of the cold gas 
shows that CO and Hi are displaced west of the center 
of NGC 4438, although CO is presen t also in the center 
of NGC 4435. iVollmer et al . ( 2005) proposed a model 
of the encounter based on the combined picture coming 
from the ISM observations and simulations. These lat- 
ter authors indicate that the evolution of NGC 4438 has 



1 INAF Osservatorio Astronomico di Padova, 
vicolo delF Osservatorio 5, 35122 Padova, Italy; 
pasquale.panuzzo@oapd.inaf.it, alessandro.bressan@oapd.inaf.it, 
lucio.buson@oapd.inaf.it, marcel.clemens@oapd.inaf.it, gian- 
luigi.granato@oapd.inaf.it, roberto.rampazzo@oapd.inaf.it 

2 Instituto Nacional de Astroffsica, Optica y Electronica, 
Apdos. Postales 51 y 216, OP. 72000 Puebla, Pue., Mexico; 
ovega@inaoep.mx, jvaldes@inaoep.mx 

3 INAF Osservatorio Astronomico di Trieste, Via Tiepolo 11, 
1-34131 Trieste, Italy; silva@oats.inaf.it 

4 Scuola Internazionale Superiore de Studi Avanzati (SISSA), 
via Beirut 4, 34014, Trieste, Italy; danese@sissa.it 



been influenced by its passage though the intra-cluster 
medium (ICM) , where ram-pressure stripped the H I in 
about 100 Myr. During this period, NGC 4435 passed 
through the disk of NGC 4438 at high velocity (« 800 
km s _1 ), but the ISM-ISM interaction of the two galax- 
ies had a less important effect on the final distribution of 
the gas than the cluster ram-pressure stripping exerted 
by the ICM. 

If the above is the general picture of the NGC 4435/38 
encounter, its effects on the evolution of NGC 4435 are 
less cl ear. In low density environments. iDominigue et alJ 
( 2003) showed that early-type galaxies show little or no 
IR or optical signatures of star formation when physically 
paired to a gas rich companion. However, some of the 
early-type members do show significant levels of IR and 
Ha emission that indicate both thermal (star formation) 
and non thermal (AGN) activity. These pairs are among 
the best candidates for direct interaction fueling of both 
starburst and active galactic nuclei. 

NGC 4435 currently has an ISM. HST images reveal 
the presence of a c ircumnuclear disk of A" radius (360 pc, 
inclin ation 44 deg: lHo et al.l2002tlCaDetti fc Balmaverdei 
2005). Moreover, n ebular emission lines were detected 
by se veral authors Kennev et alJ 119951: IHo et al.l 119951 
11997ft : in particular IHo et alJ ljl997ft . using the classical 
diagnostic diagrams, classified this galaxy as T2/H:, i.e. 
a transition object bet ween LINER 2 and H II. 

iCoccato et al.l (12006: studied the ionized gas kinemat- 
ics using HST-STIS in the central 2" looking for a super 
massive black hole (SMBH). The rotation curve is very 
symmetric, indicating a regular rotation. Their anal- 
ysis indicates that the mass of the SMBH (upper limit 
7.5 x 10 6 M Q ) is lower than that expected from the Mbh- 
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a and MsH-^buige relations. Coccato et al. suggest that 
the SMBH i n NGC 4435 cou l d be o f the "laggard" type as 
discussed in lVittorini et alJ l)2005|) . although these seem 
to lurk among galaxies that spend most of their lifetime 
in the field, at least apparently at odds with NGC 4435's 
cluster location. 

In this paper we study the Spitzer IRS MIR spectra 
of NGC 4435 presented in lBressan et al.l l)2006|) and per- 
form a comprehensive analysis of the entire SED of the 
galaxy. In particular, the study of the MIR region, rich 
in spectral features (emission lines, PAH, dust emission 
etc.) provides an independent vie w of the nature of 
the nuclear activ ity of this galaxy (Genzel et all I1998L 
iLutz et al.lll998t iRoussel et atl I2001D . Combined with 
the analysis of the panchromatic SED, it allows the char- 
acterization of the star formation activity and the esti- 
mate of the mass of stars formed during this interaction, 
shedding light on the rejuvenation mechanism in inter- 
acting early-type galaxies. 

The paper is organized as follows. In §[5] we provide 
details of the observations and data reduction methods. 
In §|3]we describe the results of the study of the spectrum 
in the MIR region, while in §0]we extend the modeling to 
the entire galaxy SED. In §[S]we discuss the result of our 
analysis and compare them with other studies. Finally 
in §0we summarize the conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. IRS spectra 

Spitzer IRS spectral observations of NGC 4435 were 
obtained during the first Spitzer General Observer Cy- 
cle on 2005 June 1 as part of a program (ID 3419, 
iBressan et al.l2006h to study early-type galaxies that be- 
lon g, to the colour - magnitude relation of the Virgo clus- 
ter llBower et al"1ll992D . 

The observations were performed in Standard Staring 
mode with low resolution (R ~ 64-128) modules SL1 
(7.5-15.3/xm), SL2 (5-7.6/xm) and LL2 (14.1-21.3/im), 
with 6 exposures of 60 seconds for each SL module and 
10 exposures of 120 seconds for the LL2 module. 

Since the adopted IRS pipeline (version SI 2) is specif- 
ically designed for point source flux extraction, we have 
devised a new procedure to flux calibrate the spectra that 
exploits the large degree of symmetry that characterizes 
the light distribution in early-type galaxi es. The reduc- 
tion p rocedure is described in detail in Bressan et al. 
(2006) thus here we only summarize the main steps. 

For SL observations, the sky background was removed 
by subtracting observations taken in different orders, 
but at the same node position. LL segments were sky- 
subtracted by differencing the two nod positions. 

We obtained new e _ /s to Jy flux conversions by ap- 
plying a correction for aperture losses (ALCF) and a 
correction for slit losses (SLCF) to the flux conver- 
sion tables prov i ded b y the Spitzer Science Center (e.g. 
iKennicutt et alJ [2003) . By applying the ALCF and 
SLCF corrections, we obtained the flux received by the 
slit. 

For each galaxy of the program, we simulated the cor- 
responding observed linear profile along the slits by con- 
volving a wavelength dependent bidimensional intrinsic 
surface brightness profile with the instrumental point 
spread function (PSF). The a dopted profile is a two di- 
mensional modified King law l|Elson et al.11198'^1 . By fit- 




FlG. 1. — HST WF PC2/F814W image of the NGC 4435 nucleus 
ICoccato et al. 2006) with IRS SL and LL slits (dashed lines) dur- 
ing our observations. The rectangles show the aperture used to 
extract the spectrum. The dot-dashed circle shows the 5" circular 
aperture used to define the SED of the central region from broad 
band data (see Sec. [1J. 



ting the observed profiles with the simulated ones, we 
can reconstruct the intrinsic profiles and the correspond- 
ing intrinsic SED. This procedure has also the advantage 
to recognize whether a particular feature is spatially ex- 
tended or not. 

Finally the spectrum was extracted in a fixed width 
(18" for SL and 10.2" for LL) around the maximum 
intensity. In order to account for the different extrac- 
tion area, the LL spectrum has been rescaled to match 
that of SL (3.6"xl8"). In Figure [T] we overplot the 
extraction area and the IRS slit positions on an HST 
WFPC2/F814W im age of the NGC 4435 nucleus by 
ICocca.toet alJEonfl. 

The uncertainty on the flux was evaluated by consid- 
ering two sources of noise: the instrumental plus back- 
ground noise and the poissonian noise of the source. 
The former was evaluated by measuring the variance of 
pixel values in background-subtracted coadded images 
far from the source. The poissonian noise of the sources 
was estimated as the square root of the ratio between 
the variance of the number of e~ extracted per pixel in 
each exposure, and the number of the exposures. The 
total noise was obtained by summing the two sources in 
quadrature and by multiplying by the square root of the 
extraction width in pixels. We notice that the overall 
absolute photometric uncertainty of IRS is 10%, while 
the slope deviation within a single segment (affecting all 
spectra in the same way) is less than 3% (see the Spitzer 
Observer Manual). 

2.2. IRAC and MIPS images 

IRAC and MIPS images of NGC 4435 were retrieved 
from the SSC archive. They were taken as part of the 
program ID 3649. Fluxes at 3.6, 4.5, 5.8 and 8 /im 
were measured from the IRAC images as output from 
the "post-BCD" pipeline. The proximity of NGC 4435 
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to NGC 4438 results in a weak contaminating flux com- 
ponent. This was treated by estimating the median back- 
ground flux from the image containing NGC 4435 so 
that the subtracted background included the contam- 
inating light. Fluxes were extracted in circular aper- 
tures centered on the brightest pixel at each wavelength. 
Given that the calibration of extended sources for IRAC 
is still somewhat uncertain the flu xes so extracted have 
errors of - 10% (|Dale et al.112005^ . Similarly, we mea- 
sured fluxes at 24, 70 and 160 /jm from MIPS mosaics 
as output from the "post-BCD" pipeline. Fluxes were 
extracted in circular apertures centered on the brightest 
pixel at each wavelength, subtracting the median back- 
ground flux measured far from the sour ce. Uncertainties 
for MIPS fluxes are quoted to be ~ 20% i|Dale et al.l2005l 
see also the Spitzer Observer Manual). It is worth not- 
ing that in all MIPS images the galaxy appares almost 
as a point source; in particular at 24 /jm where the angu- 
lar resolution is higest the galaxy image is only slightly 
larger than the instrumental PSF (6"). 

3. ANALYSIS OF THE MID-IR SPECTRUM 

The final flux calibrated IRS spectrum of NGC 4435 
is shown in the top-left panel of Figure |31 The spec- 
trum is dominated by numerous PAH features and shows 
several nebular emission lines, [Arli]7/tm, [Neil]12.8/im, 
[Ne in]15.5/im and [S iii]18.7/im, and molecular hydrogen 
rotational lines, H 2 S(3) 9.66/tm, H 2 S(l) 17.04/mi. 

In order to measure the fluxes of emission lines and 
to identify the PAH features, we fitted the spectrum via 
decomposition into various components: PAH emission 
features, emission lines and underlying continuum. 

For all the PAH features but those at 11.2 and 12.5 
/im, we assumed a Drude profile. The features at 11.2 
and 12.5 /im are ve ry asymmetr i c thus we have used the 
profiles observed by Ho nv et all (|200lD . convolved to the 
SL1 resolution. I t is worth noticing that the emission 
profiles derived bv lHonv et alJ (|2001) do not account for 
the blending of the wings of the PAH features. Conse- 
quently our underlying continuum shows a broad feature. 
However, since the main purpose of this section is to es- 
timate the emission line intensities, this does not affect 
our conclusions. 

For the emission lines we adopted a Gaussian profile, 
with a FWHM given by the resolving power of the corre- 
sponding IRS module (see the Spitzer Observer's Man- 
ual, v. 6.0, sec. 7.1.6). The resulting fits of the three 
most interesting spectral regions are shown in Figure |3J 

3.1. PAH features 

In the IRS spectrum of NGC 4435 we can identify many 
PAH emission features already known from ISO and pre- 
ISO studies: at 6.2, 7.7 8.2, 8.6, 11.2, 11 9, 12.7, 13.5 
14.5, 15.9, 16.4 /mi (e.g. iLi fe Dramll200lT: ISturm et all 
2000). Besides those features, we also see a number 
of new components that have been revealed by recent 
Spitzer observations. These are the feature at 10.7 /im 
and, most notably, the complex around 17 /an. The na- 
ture of the fo rmer is unknown but there is similar evi- 
dence in M82 jForster Schreiber et al.ll2004D. NGC 7714 
l|Brandl et aTll2004[) and NGC 7331 l|Smith et ai]l2004]) 
spectra. The complex around 17 /im (Figure bottom- 
right panel) is similar to that observed in the spiral 
galaxies NGC 7331 and in the extended emission of M82 
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Emission line intensities. 
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tensity provided by the fitting procedure, how- 
ever its value is unreliable because it falls in 
the wing of the broad PAH 7.7 /im emission 
feature. 

ijEneelbracht et al] [2006) . This complex is superimposed 
on H 2 S(l) emission but in our case, since the H2 S(l) 
emission is fainter than that seen in NGC 7331 and M82, 
we can more accurately ext ract the underlying emission 
components. In fact, while ISmith et aD l)2004|) claimed 
the presence of a new broad emission feature at 17.1 /im 
with a FWHM of 0.96 /im, in our spectrum this fea- 
ture appears to be a blend of two components; a nar- 
row component centered at 17.4 /tm with a FWHM of 
0.33 /tm and a broad component centered at 16.9 /tm with 
a FWHM of 1.1 /tm. Notice th at the component cen- 
tered at 17.4 /im was observed bv lVan Kerckhoven et alJ 
(2000) in various galactic star formation regions. Thus 
our observations indicate that the new PAH emission 
feature is centered at 16.9 /im. In the same spectral 
region we also find an emission feature at 19.07 /tm 
th at was observ e d in the NCG 7023 reflection nebula 
bv lWerner et alJ l)2004|) . 

3.2. Emission lines 

In the MIR spectrum of NGC 4435 we detect only low 
ionization emission lines. The corresponding intensities 
are listed in table ^ The reported uncertainties in the 
line fluxes are la as derived from the line fitting. For 
some interesting higher excitation emission lines we can 
only provide 3<r upper limits. In the following we discuss 
the main implications of our detections. 

3.2.1. Constraints on AGN excitation 

Nebular emission lines provide a powerful diagnos- 
tic tool to determine the nature of the source of en- 
ergy powering the emission. Classically, the determina- 
tion is done using optical emission lines. Optical spec- 
tra of_thenuclear region of NGC 4435 were observed 
bv iHo e t al.l (IT 995I). Using t he diagno s tic di agrams of 
IVeilleux fc Osterbrockl l(T987|) . iHoltall lfT997|) classified 
NGC 4435 as a transition object between LINER and 
H11, via the value of [Oi]6300/Ha (=0.13) which, how- 
ever, has an uncertainty of 30-50%. Revisiting the 
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Fig. 2. — Spitzer IRS MIR spectrum of NGC 4435. Top-left: entire observed spectrum; arrows indicate the positions of PAH features, 
nebular emission lines and H2 rotational lines. Top-right: the 6.5— 10/xm spectral region split in the different components. The error 
bars correspond to 1 a noise fluctuations. The lower smooth solid line corresponds to the continuum component; the dotted line to the 
continuum plus PAH features, the dot-dashed lines show the continuum plus emission lines. Finally the upper solid line represents the sum 
of all components. Lines for which we could provide only upper limits are included in parenthesis. Bottom-left: similar to top-right panel 
but for the 10-14. 6^m spectral region. Bottom-right: similar to top-right panel but for the 15-20/^m spectral region. 



cla ssification using the AGN/H11 separation proposed 
bv IKewlev et al.l l)2001|) . the [Oi]6300/Ha ratio is only 
marginally larger than that expected for H11 regions. 

The IRS spectrum allows us to make use of the ratio 
[Ne in]15.5/[Neii]12.8 which is a strong infrared diagnos- 
tic for the hardness of the ionizing source. We compared 
the observed ratio of [Ne Hi] 15.5/[Nei i1 12.8=0.38 with 
CLOUDY models ijFerland et al.lll99c| . We computed 
several photoionization models using as ionizing source 
the spectrum of a 35000 K star (with Kurucz atmo- 
sphere model) plus a 10 K blackbody, repre senting the 
AGN contribution (a similar test was done bv lLutz et alJ 
1998). Assuming an ionization parameter log(J7) = —2, 
and nn = 70 (derived from Sll optical lines), the ob- 
served ratio [Ne in] 15.5/ [Neil] 12.8 limits the AGN con- 
tribution to the ionizing flux to be less than 2%. 

3.2.2. H2 rotational emission 

In Table fj] we provide the intensity of detected H2 ro- 
tational emission lines. 

The first excited rotational state of H2 that can provide 
a transition toward the fundamental level has an angular 
moment J — 2 and an energy of 510 K, thus the H2 does 
not radiate for temperatures typical of the cold molecular 
medium (some 10 K) but only for temperatures above 
100 K. The investigation of the state of warm molecular 
gas can be done using the excitation diagram, shown in 
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Fig. 3. — Excitation diagram for H2 rotational lines observed in 
the NGC 4435 IRS spectrum. Diamonds show the natural loga- 
rithm of the number of molecules in the state i computed with eq. 
fl] divided by the statistical weight plotted versus the energy of the 
excited state with respect to the ground state. Error bars represent 
ler errors on line intensity, while arrows show upper limits on H2 
S(4) and H2 S(2) lines. The dashed and dotted lines show the fits 
to the warm and hot components, respectively. 



Figure I?. 2. 21 where the natural logarithm of the number 
of molecules in the state i (iVj) divided by the statistical 
weight of the state (<?i) is plotted against the energy (Ei) 
of the excited state with respect to the ground state. The 
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number of molecules in state i is given by (assuming a 
density smaller than the critical density): 



abundance relative to hydrogen is given by: 



Ni = 



(1) 



where Li is the luminosity of the line emitted by the 
transition starting from the state i, Ai is the transition 
probability, and Vi the frequency of the transition. 

Assuming local thermal equilibrium, the ratio between 
number of molecules in state i and the total number of H2 
molecules per unit of volume is given by the Boltzmann 
equation (jParmar et al.lll99"l"|) : 



7lH, 



hcB 



2kT 6XP \ kT 



E, 



(2) 



where B is the H2 rotational constant. 

For an isothermal gas all the data in the excitation 
diagram of Figure 13.2.21 should lay on a linear relation 
y = a + bx with b = — 1/T. However, the observed val- 
ues do not lay on the same line, suggesting the presence 
of two (or more) components at different temperatures, 
a "warm" component at T < 240 K traced by the S(l) 
line (and the S(2) upper limit) and a "hot" component 
at T ~ 630 K, traced by lines at shorter wavelengths. 
With equations ^] and [21 we can estimate the mass of the 
"warm" and "hot" components. Assuming a tempera- 
ture of 200 K for the warm component, we find a mass 
of 6.3 x 10 4 M Q , while for the "hot" component we find 
a mass of 7.1 x 10 3 M Q . 

The excitation mechanism for hot and warm molecular 
gas in galaxies is commonly attributed to UV flux from 
young stars in photodissociation reg ions (PDR), shocks 
or X-ray flux from an a ctive nucleus (|Burton et alJll992T 
IRigopoulou et al.H2002|) . 

The upper limit on the temperature derived from 
S(l) and S(2) lines exc ludes AGN excitation. In fact 
IRigopoulou et af] ()2002|) find that in Seyfert galaxies the 
temperature from these lines is typically around 350 K. 
Indeed, the temperature we find is compatible with PDR 
models by iBurton et alJ ©92). Finally, it is worth 
noting that the mass of the warm gas is only a tiny 
fraction (~ 0.06%) of the total molecula r gas observed 
in NG C 4435 from CO lines (~ 1O 8 M , iVollmer et all 
2005). This fraction is also significantly smaller than 
that found in typical powerful starbursts (1 - 10%) o r 
in Seyfert galaxies (up to 30%. IRigopoulou et alJ l20021. 
This is consistent with our finding that the star formation 
in NGC 4435 has already declined and is more charac- 
teristic of a post-starburst phase, as we discuss in the 
following sections. 

Note that the above values are not corrected for dust 
attenuation. However, the attenuation we derived from 
the model fit to the global SED in the following sections 
is quite low (A\ ~ 0.1 — 0.2) at these wavelengths. 

3.2.3. Abundance determinations 

From the intensity of the Pfa, neon and argon lines 
we could obtain an estimate of the metallici ty of the ion- 
ized g as, following the procedure outlined in lVerma et al.l 
(2003). In fact, for a generic element X and as long as 
the current density is lower than the critical density for 
collisional de-excitation of the involved transitions, the 
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where jpf a is the emissivity for the Pfa, and <7[xii] is 
the collisional excitation rate of [X 11] transition (for elec- 
tronic density near to zero). The above abundances can 
be translated into metallicity estimations by dividing 
them by standard ISM abundances. 

We notice however that, while we are confident on the 
determination of metal line intensities, the Pfa falls in 
the wing of a broad PAH emission feature. The mea- 
surement of Pfa intensity is thus unreliable, even if the 
fitting procedure gives a small uncertainty (reported in 
table nj. For this reason we do not proceed with a direct 
determination of the abundances in this way. However 
an alternative way to exploit the good measurements of 
metallic lines in order to derive the metallicity is to use 
the hydrogen line intensity as predicted from the fit to 
the global SED of the galaxy. This will be discussed in 
a following section. 

4. ANALYSIS OF THE PANCHROMATIC SED 

The analysis of the MIR spectrum made in the pre- 
vious section indicates that a star formation activity is 
present in the circumnuclear disk of NGC 4435. In order 
to better characterize the emission mechanisms and the 
stellar population content of this galaxy we will model 
the panchromatic SED of NGC 4435 by means of our 
spectro-photometric code GRASIL 5 with the recently 
updated treatment of PAH emission based on ISO data 
(for details on GRA SIL seelSilva et al.lll99St ISilvalll999t 
iGranato et al J 1200(1 iBressan et aTT T2002: iPanuzzo et alJ 
l2003ll2005UVega et al.H2005j) . 

The panchromatic SED was constructed by combin- 
ing the IRS spectra (including the main nebular lines) 
with archival data matching as far as possible the IRS 
apertures. J, H and Ks fluxes within an aperture of 5" 
radius were taken from the 2MASS database at IRSA 6 . 
IRAC fluxes in the 3.6, 4.5, 5.8, and 8.0 /im bands were 
extracted in the same aperture from archival images as 
explained in section 12.21 IRAS 60, and 100 /im fluxes 
(obtained from NED) were used to define the FIR SED 
together with MIPS total fluxes. Although the aper- 
tures of MIPS and IRAS are larger than the IRS aper- 
ture, the above fluxes are dominated by emission from 
the dusty circumnuclear disk as clearly shown by the 
24 /im MIPS image. This is not true, however, at 12 
/im, where we expect a significant contribution from 
the old stellar population; thus the 12 /im IRAS flux 
has not been included in the SED. We also excluded 
the IRAS 25 /im flux (0.22 Jy, from NED) because it 
is not compatible with the 24 /im MIPS flux; this is 
probably due to an incorrect background subtraction 
due to the poor resoluti o n of I R AS. Radio data were 
taken f rom INagar et al.1 (I2000D . IWrobel fc Heenschenl 
(EMI), iBecker et all (|1991|) . iDressel fc Condon! (|1978|) . 

5 Available atlhttpV / adIibitum.oatsmiaLi t / silva/ default ,html| or 
http: / /web. oapd.mat.it/granato/grasii. html 

D The NASA/IPAC Infrared Science Archive is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space Administra- 
tion. 
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TABLE 2 

Broad band fluxes of NGC 4435 applicable 

FOR 5" CENTRAL REGION 



A 

(H 


Band / 
Instrument 


Flux 

(Jy) 


Uncert. 
(mjy) 


Ref. 


1.25 


j 


0.11 


12.1 


\ 


1.65 


H 


0.14 


15.3 


\ 


2.17 


Ks 


0.12 


13.2 


I 


3.6 


IRAC 


0.056 


5 6 


> 


1.5 


IRAC 


0.033 


3.3 


2 


5.8 


IRAC 


0.039 


3.9 


2 


8.0 


IRAC 


0.062 


6.2 


2 


23.7 


MIPS 


0.11 


22.0 


2 


60 


IRAS 


1.99 


597.0 


3 


71.0 


MIPS 


2.27 


454.0 


2 


100 


IRAS 


4.68 


1404.0 


3 


156 


MIPS 


3.7 


740.0 


2 


20000 


VLA 


<0.0011 




4 


60000 


VLA 


>0.0012 




5 


61856 




<0.007 




6 


126050 




<0.004 




7 


200000 


VLA 


0.0022 


0.216 


8 



References. — (1) IRSA; (2) this work; (3) 
NED; (4) Nagar et al. 2000; (5) Wrobel & 
Heeschen 1991; (6) Becker et al. 1991; (7) Dres- 
sel& Condon 1978; (8) FIRST VLA survey. 



and the FIRST- VLA survey. The SED is completed 
with GALEX UV fluxes extracted within a 5" aperture 
(Boselli & Cortese, private communication). The broad 
band fluxes are given in Table |2j 

4.1. The models 

Since the broad band SED obtained in the previous 
section is produced by the star-forming circumnuclear 
disk and the old population, including also the stars seen 
in projection, we adopted a simple composite model for 
the fit: a starburst plus an old stellar population. The 
old stars are modelled as an unattenuated simple stellar 
population (SSP) with given age and metallicity. The 
SSPs are based on Padova stellar models and include the 
contribution from dusty circumstellar envelopes around 
asymptotic giant branch stars. This fact may be relevant 
because these old stars leave a n excess around 10 /jm 
ijBressan et al.lll998l l200lU2006j) . 

The starburst contribution was modelled with 
GRASIL assuming different values for the age (it>), 
metallicity and star formation history. The latter was 
assumed as an exponentially decreasing SFR parameter- 
ized by the e-folding time, tv We consider two phases 
for the ISM: the molecular clouds (MCs) which surround 
young star clusters and the diffuse component. For the 
distribution of both the diffuse ISM and stars we adopted 
a disk geometry, with an inclination of 45° with respect 
to the line-of-sight (as suggested by observations, e.g. 
iCoccato et alJ l2006j) . The density of dust and stars is 
assumed to be exponential in both radial and vertical 
direction, with different scales for dust and stars (see 
ISilva et alJll99l for more details). 

Among the GRASIL physical quantities that mostly 
affect the shape of the SED we recall (i) the total gas 
mass M gas , (it) the fraction of molecular mass 7 to total 
gas mass (M mo i/M gas ), (Hi) the optical depth of molecu- 

7 In the GRASIL model we identify the molecular gas with the 



TABLE 3 

Parameters for the burst model library. 





log(tb) 


tcsc A^mol/A^gaa T\ 


(Myr) 


(y) 


(Myr) 


10 - 100 


6.0 - 8.8 


1 - 99 0.10 - 0.99 0.01 - 180 



TABLE 4 

Best fit parameters for the burst model. 







z t csc Af mol /A/ gas n 


(Myr) 


(Myr) 


(Myr) 


55 


186 


0.02 3.0 0.20 0.25 



lar clouds at 1/xm (ti), and (iv) the escape time-scale of 
newly born stars from their parent MCs (t esc ). The dust 
composition and optical properties are set to reproduce 
the extinction and emission properties observed in the 
diffuse ISM of our Galaxy. The gas-to-dust mass ratio 
(G/D) is set proportional to the metallicity (Z), with 
G/D = 110 for solar metallicity. 

For the nebular li nes we adopted t he lar ge Hii region 
library computed bv iPanuzzo et al. (2003). Finally, the 
adopted IMF is a Salpeter IMF l)Salpeterlll955|) in the 
mass range from 0.1 to 100 M@. 

4.2. Results of the SED fit 

We computed a large grid of simulated SED with 
GRASIL. In Table we list the ranges of values used 
for the library. Among these models, we identified the 
best fitting model as the one that maximize a likeli- 
hood function. The likelihood function was built using 
UV, MIPS, IRAS and radio data and their uncertainties 
quoted in Table 2 and th e intesities of observed line Ha, 
H/3 (from lHo et aT]|l997|) and [Neil]12.8/xm. The best fit 
model is shown in Figure |4~2*1 where the left panel shows 
the panchromatic fit, from the UV to radio wavelengths, 
while the right panel is a zoom into the MIR spectral 
region. The corresponding parameters and some impor- 
tant derived quantities are displayed in Tables 0] and [5] 
respectively. 

The wide wavelength coverage of the SED imposes 
strong constraints to the model, because different spec- 
tral regions probe the age sequence of the stellar popula- 
tions. Broadly speaking, nebular emission lines are dom- 
inated by massive stars, radio emission is contributed by 
stars younger than about 50 Myr while the FIR SED in- 
cludes also a significant contribution from intermediate 
mass stars. The UV/FIR ratio implies that the burst 
luminosity is almost completely reprocessed by dust and 
emitted in the FIR, thus the FIR luminosity fixes the 
bolometric luminosity of the burst. The emission line lu- 
minosities and radio fluxes constrain the star formation 
rate in the last 50 Myr determining the timescale t\,. 
Once the current SFR and Th are determined, the age of 
the burst is fixed by its luminosity. Finally, the Ha and 
H/3 constrain the attenuation at optical wavelenght (see 

clouds surrounding young star clusters. Molecular gas could also be 
present in the diffuse component but the code doesn't differentiate 
between atomic and molecular gas in this phase. 
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TABLE 5 

Derived quantities of the SED models. 



(Mq) 


(Mq) 


(Mq) 


L SB 

(96) 


L o\d 

(96) 


(SFR) 
(M yr- 1 ) 


SFR C 

(Mq yr- 1 ) 


/I mod 
Ay 

(mag) 


1.40 x 10 s 


1.36 x 10 s 


8.2 x 10 9 


32.0 


68.0 


0.75 


0.089 


2.69 




1 2 3 

log(X) [Mm] 




-2.0 



10 15 
X [/tim] 



Fig. 4. — Comparison between the observed SED of the central region of NGC 4435 and our model. The thick solid black line represents 
the model for the total SED, i.e. the starburst component plus the old stellar component; the three dots-dashed darkgreen line represents 
the contribution from the old stellar population, and the dashed blue line represents the total contribution from the burst of star formation, 
the dotted red line represents the emission from molecular clouds. The dot-dashed green line represents the diffuse medium emission and 
the thin cyan solid line denotes the emission from stars of the starburst component without applying the extinction from dust. The red- filled 
circles are the broad band data. Left: Comparison from 0.1/im to 100 MHz. Right: Comparison for the MIR wavelengths. The thickest 
blue solid line represents the IRS Spitzer spectrum. 



section l5~2*l for more details). 

Our assumption of an exponentially declining SFR is 
the simplest model consistent with a Schmidt-type law 
that reproduces the multiwavelength SED. Other more 
complex scenarios could probably be able to fit the data 
introducing, however, further free parameters that are 
not justified. Indeed, stars in the obscured burst older 
than ~ 50 Myr do not provide observable "features" to 
disentangle the detailed shape of the star formation his- 
tory. 

It is also interesting to consider the effects of a underly- 
ing star formation activity already present in NGC 4435 
before the interaction (i.e. older than ~ 200 Myr). Stars 
originated from this activity would unlikely remain con- 
fined to the dusty disk. For example in our galaxy, stars 
older than 150 Myr show larger spatial distribution tha n 
younger stars (R obin et alJl2003t iPanuzzo et alJ feOOG). 
Thus the intermediate age stars would now be mostly 
outside the dusty disk and, as a consequence, they would 
suffer from only a small or negligiable attenuation. Since 
stars 200-500 Myr old are still important emitters at UV 
wavelengths, they would contribute significantly to the 
GALEX fluxes. However, the measured UV GALEX and 
NIR fluxes within the 5" aperture (Fig. 14.21 left panel), 
are already entirely accounted for by the old population 
and the obscured starburst. Actually, the observed FUV 
flux within the 5" aperture corresponds to an activity not 
larger than SFR - 3 x 10" 4 Af Q yr" 1 between 200 and 
500 Myr ago. This suggests that no significant under- 
lying intrinsic star formation activity was present and, 
consequently, that the observed burst was triggered by 
the interaction. 



The old stellar population dominates the emission in 
the optical and NIR range but is significant also in the 
MIR. In fact we see from Figure l4~2l that the fractional 
contribution of the starburst to the IRS spectrum in- 
creases with wavelength. This is consistent with the re- 
sults of the spatial analysis along the slit presented in sec- 
tion l2.ll We find that the continuum below ~ 6 fim looks 
extended, while the emission at longer wavelengths (in 
particular the PAH features) looks unresolved, in agree- 
ment with the different concentration obse rved in ISO 
LW2 and LW3 images (|Boselli et al.l 1200^) . IRAC and 
MIPS images, as detailed below. 

Finally we notice that the fit of the SED does not re- 
quire an AGN component, which adds further support 
for the lack of AGN activity in NGC 4435. 

In the following we discuss in more detail the results 
of the fit. 

5. DISCUSSION 

5.1. The star formation history and starburst 
luminosity 

The old stellar component is well fitted by a solar 
metallicity SSP of 8 Gyr and a mass of 8.2 x 10 9 M Q . 
Note that the value here deri ved for the stel l ar ma ss de- 
pends on the adopted IMF. iCoccato' eta!] p006|) give 
an estimate of the mass to light ratio (M/L) of the 
population that can be directly compared with our old 
stellar population. Using the / band they measured 
M/L(I) = 2.2-M Q /L @ (I), while for our old stellar popu- 
lation it is M/L(I) = 2.98 -M & /Lq(I). The discrepancy 
in the M/L can be solved by changing the slop e of the 
IMF at lower masses (as in the Kroupa IMF, iKroupal 
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1998) or adopting a higher lower mass limit (e.g. 0.2 
Mq) without changing the predicted SEDs. 

The nuclear disk is in a post-starburst phase (jh = 55 
Myr and % — 186 Myr) with a current star forma- 
tion rate (SFR C = 0.089 M Q yr _1 ) about one order 
of magnitude smaller than the average star formation 
rate ((SFR) = 0.75 M Q yr" 1 ). The quality of radio 
data, however, prevents a determination of Tb better than 
±20%; this results in an uncertainty on the age of the 
sturburst t^ of the same order. The mass of gas converted 
into stars during the burst is MfJ 5 = 1.40 x 10 8 M ; 
however with the adopted IMF, only ~ 80% still survive. 
The star forming component emits 32% of the bolometric 
luminosity in the central 5" of the galaxy, while account- 
ing for 1.4% of the stellar mass. Most of the starburst 
luminosity is provided by A stars, reprocessed by the dif- 
fuse dust and re-emitted in the 8-1000 /im spectral region 
(L m ~ 5.5 x 10 35 W). 

A population of A stars in the central ~ 550 
pc has already bee n detected by o ptical obser v ations 
(jKennev etaflll995|) . More recently, ISarzi et alJ (|2005f) 
estimated that the contribution of the young burst from 
3000 A to 5700 A in the nuclear region (0'/13 radius) is 
around 12%. In the same spectral region we estimate 
that the fractional contribution of the starburst is only 
6%. It is however hard to compare the two estimates 
because our value refers to a much larger region (5"). In- 
deed it is sufficient that in the very cen tral region the at- 
tenuation be reduced (as estimated bv ISarzi et alJl200"5^ 
to eliminate this discrepancy without affecting signifi- 
cantly the constraints from the other spectral bands. 

The strongly declining SFR makes this galaxy a typ- 
ical cases where different SFR indicators, probing dif- 
ferent timescales, provide estimates that are apparently 
discrepant. In fact with the classica l SFR estimation 
based on the IR luminosity provided bv lKennic utt (1998) 
we get SFRir = 0.25 Mq yr _1 , in between the av- 
erage SFR a nd the current one. Using the 1.4 GHz 
flux with the Pa nuzzo et alJ l)2003|) calibration we get 
SFRi.4GHz — 0.117 Mq yr _1 . This value is more similar 
to the current SFR value than the IR based estimate be- 
cause the radio emission is due to free-free emission from 
HII regions and synchrotron emission from core collapse 
supernovae. 

5.2. The geometry and the attenuation law 

Though the star formation history is well constrained 
by the fit of the total SED, a further refinement to the 
relative spatial distribution of stars and dust can be ob- 
tained by the comparison of the predicted and observed 
Ha and H/3 line intensities. It was shown in previous 
works that a system with an equal spatial distribution 
of stars and dust produces a much smaller reddenin g 
than a "screen" geometry (e.g. iWitt fc Gordonl [2000). 
With the assumption of equal scales for stars and dust 
and given the strong attenuation implied by the ob- 
served UV/FIR ratio, the Ha intensity predicted by our 
best fit model would be much small er than the obse rved 
value (7(Ha)= 2.82 xl0~ 17 W m- 2 . IHo et aT!IT997h and 
the predict ed gas mass about twice the quoted value 
~ 10 8 M Q l|Vollmer et all 120051) . To solve this discrep- 
ancy we need to decrease the attenuation in the optical 
band leaving unchanged that in the UV. Without mod- 
ifying the optical properties of dust grains, this can be 
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Fig. 5. — Comparison of Cardelli extinction curves with different 
Rv with the attenuation curve derived from the best fit model 
(solid line). 



achieved by assuming a geometry where the diffuse dust 
is more extended than the stellar distribution. This has 
the effect of increasing the reddening, because it reddens 
the attenuation law making it more similar to the extinc- 
tion of a uniform slab. Adopting a radial (vertical) scale 
length of 160 (22) pc for the diffuse dust and of 70 (14) 
pc for the stellar component the predicted intensities for 
Ha and H/3 are I(Ha)= 2.70x 10" 17 W m" 2 , and /(H/3)= 
4.99xl0 -18 W m -2 well in agreement w ith the observe d 
value (4.99xl0" 18 W m" 2 for H/3 from iHo et aD 11997ft . 
The required dust mass is 1.23 x 10 6 M Q . For a solar 
metallicity, as predicted by other observables discussed 
below, this implies a gas mass Mq ~ 1.36 x 10 8 Mq, 
fairly near to the observational value. 

The estimated mean attenuation of the system in the 
V band is ^4y od = 2.69 mags. The resulting attenuation 
law is shown in Figure [3] In the same figure we plot, 
for comparison, thre e different extinction curves from 
iCardelli et al.1 l)1989j) . We stress that the dust proper- 
ties adopted in the model are such that they reproduce 
the Galactic extinction law, corresponding to i?v=3.1. 
Thus the geometry of the model (including the effects 
of age dependent attenuation) produces an effective at- 
tenuation law which is very different fro m the adopted 
intrinsic extinction properties of the dust (jGranato et alJ 
1200(1 IPanuzzo et al]l2(Ml200fil) . 

5.3. The gas and star metallicity 

The best fit to the NGC 4435 SED is obtained by 
adopting solar metallicity for both stars and gas in the 
starburst. 

Foll owing the detai l ed stu dy of nebular emission mod- 
els of IPanuzzo et al.l l|2003|) we can also compute the 
predicted intensities of selected emission lines. This is 
particularly useful for abundance determinations with 
infrared lines in the absence of g ood measurements o f 
hydrogen recombination lines (see IPanuzzo et aLl l2003) . 
as in the present case. In fact from the SFR obtained 
from the fit of the SED, GRASIL also computes the ion- 
izing flux and thus, the observed intensity of the main 
emission lines. This can be considered as an alternative 
method to obtain the ionizing flux from a direct mea- 
surement of the hydrogen recombination lines. Adopting 
solar metallicity and standard ISM abundances we obtain 
the following predicted intensity: J(Pfa) = 3.46 x 10~ 18 
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W m -2 , /([Neil]) = 7.01 x 1CT 17 W m" 2 , /([Nelll]) = 
2.60xlCT 17 Win" 2 , /([Am]) = 3.17x 10~ 17 W m~ 2 and 
/([Sll]) = 1.71 x KT 17 W m~ 2 . For the metallic lines 
these figures are in good agreement with the observed 
values reported in Table ^ thus providing an indepen- 
dent confirmation that the gas metallicity is solar. The 
value found for Pfa is much lower than that derived from 
the fit of the line, suggesting that the blend with PAH 
prevents an accurate measurement (see section 13. 2. 3[1 . 

5.4. PAHs and IRS continuum 

The left panel of Figure 14.21 shows a detailed compar- 
ison MIR Spitzer spectrum and the model. The contin- 
uum and PAH features up to 13 /an are well fitted by 
the starburst model. At wavelengths longward of 13 /im, 
while the predicted continuum matches fairly well that 
observed, the predicted PAH features are in disagree- 
ment with those observed, both in central wavelength 
and strength. Only the 16.4 /an feature is reproduced. 

We recall that our m odel for PAH emission is based on 
the ILi fc Dramel 1)200 lj) model, which predicts PAH fea- 
tures at 3.3, 6.2, 7.7, 8.6, 11.3, 11.9, 12.7, 16.4, 18.3, 21.2, 
and 23.1 /im. The model of PAH features at A < 13 /im 
was based on ISO data for the diffuse galactic medium. 
While the FWHM and the cross sections for the features 
at A > 13 /im were taken f rom the theoretical st udy of 40 
families of PAHs driven bv lMoutou et alJ l)1996|) . because 
of the lack of corresponding observational data. 

More recently, ISO extended the observed wavelength 
range, and confirmed the feature at 16.4 /jm, and also 
found new features at 14.2 and 17 .4 ^tm, as well as a 
weak plateau betw e en 15 and 20 /zm (IMoutqu et al.l2000t 
ISturm et alJ l2000l iVan Kerckhoven et all l2000|) . This 
plateau is attributed to a blend of many emission fea- 
tures provided by the interstellar PAH families present 
in the sources, and its variation can be accommodated b y 
variations in the PAH populations ijPeeters et al.l20 04a'). 
Lately, the higher sensibility of Spitzer confirmed the ex- 
istence of PAH features at 14.2, 16.4 and 14.7 /im, and 
a new and promi nent feature at 17.1 /zm, among other 
weaker ones fe.g. iBrandl et alJl2004l iSmith et"aLll2004t 
IWerner et al.l2004jr However, the lack of features at 18.3 
/im, 21.2 /mi, and 23.1 /zm has been confirmed. 

The dis agreement betw een model and data was already 
noticed bylPeeters et al.l (|2004blh ISmith et al"](|2004|) and 
IDale et alJ 1)2005(1 who, on the other hand, found very 
good agreement between the Li & Draine model and the 
data for A < 13 /im. They attributed the disagreement 
at wavelengths longward of 13 /im to the large PAH sen- 
sitivity to their global molecular structure, since they in- 
volve the motion of the molecules as a whole, and there- 
fore, they are very dependent on the exact PAH species 
in the emitting source. 

Work is in progress to improve our PAHs emission 
model based on new Spitzer observations. 

5.5. Total aperture SED 

We have also compared our models with the total aper- 
ture SED of the galaxy. The total aperture SED includes 
the U, B, V, J, H, and K bands from GOLDMINE 8 ; 
the LW2 (6.75 /im) and LW3 (15 /im) ISO bands from 

8 Galaxy On Lin e Database Milano Network, 
http:/ /goldmine. mib.infn.it 



TABLE 6 

Total aperture broad band fluxes of 
NGC 4435 



A 
(M m ) 


Band / 
Instrument 


Flux 

(Jy) 


Uncert. 
(mjy) 


Ref. 


0.2000 


UV 


3.4E-4 


0.1 


1 


0.3601 


U 


0.0219 


4.38 


2 


0.4399 


B 


0.0870 


13.05 


2 


0.5500 


V 


0.166 


24.9 


2 


1.25 


J 


0.61 


122.0 


2 


1.65 


H 


0.70 


105.0 


2 


2.17 


K 


0.53 


106.0 


2 


3.6 


IRAC 


0.32 


32.0 


3 


1.5 


IRAC 


0.18 


18.0 


3 


5.8 


IRAC 


0.16 


16.0 


3 


6.75 


LW2 


0.12 


4.61 


i 


8.0 


IRAC 


0.16 


16.0 


3 


12 


IRAS 


0.130 


29.2 


5 


15 


LW3 


0.0765 


8.73 


4 


References. - 


(1) Deharveng et al. 


2002; 



(2) Goldmine; (3) this work; (4) Boselli et al. 
2003; (5) NED. 

Note. — Complementary IRAS and radio 
fluxes were already reported in Table |21 



IBoselli et all ((2003); the UV (2000A) from Deharveng et 
al. (2002), all IRAS bands (with the exception of 25 /im), 
and the radio measurement described in section^ These 
data are complemented by our IRAC flux measurements 
and GALEX data for the full aperture. The data are 
reported in table [S] 

The fit to the total aperture SED was done simply by 
changing the amount of the old stellar population and 
using the same starbust component that fits the central 
SED, as derived in the previous sections. 

The resulting total mass of the old stellar population 
is Mj°f = 4.66 x 10 10 M . Moreover, the star forming 
component emits 8% of the total bolomctric luminosity of 
the galaxy, while its stellar mass is only 0.3% of the total 
stellar mass of the galaxy. It is worth noting that the 
starburst in N4438 (the late type companion) accounted 
for a smaller fraction (a bout 0.1%) of its total stellar 
mass (B oselli et al.ll2005l) . 

The comparison between the observed full aperture 
SED and the model is shown in the Figure El The full 
aperture SED is well reproduced at all wavelengths. We 
stress that the old stellar component provides a signifi- 
cant contribution to GALEX, ISO LW2, LW3 and IRAS 
12/im broad band fluxes. 

5.6. X-rays 

As a final step in the SED analysis we briefly consider 
the X- ray emission of NGC 4 435, observed with Chan- 
dra bv lMachacek et all (|2004l) . The X-ray spectrum was 
fitted with a thermal component in the soft (0.3-2 keV) 
band attributed to hot gas, and a non-thermal compo- 
nent responsible for the 2-10 keV band emission. The 
hard component, with a luminosity of Lx = 3.9 x 10 32 
W, is mostly concentrated within 3" from the galaxy 
center. This emission was attribu ted to a possibl e low 
luminosity AGN (as suggested by IHo et all I1997D plus 
recent star formation activity. 

The observed hard X-ray luminosity of NGC 4435 is 
about five times larger than that expected from the 
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FIR luminosity and the X-ray - FIR corre lation found 
for starbursts galaxies llRanalli et alJ 120031). Mor eover, 
adopting the calibration by Ranal li et alJ {2003), the 
SFR implied by the observed hard X-ray luminosity turns 
out SFRx-ray = 0.78 M Q yr^ 1 , one order of magnitude 
higher than the current SFR derived from the fit to the 
global SED. 

Though this seems to support the AGN nature of the 
hard emission component we need to consider the pecu- 
liar star formation history of NGC 4435. In fact, due 
to the much higher SFR at the beginning of the burst, 
the average SFR is consistent with the X-Ray estimate. 
For this reason we expect a significant contribution to 
the X-Ray luminosity from intermediate mass stars ac- 
creting onto neutron stars and/or black holes, through 
Roche-lobe mass transfer. To check this possibility we 
have computed the X-ray emission expected from the 
star formation activity derived pr eviously, with the syn- 
thesis model of iSilva et alJ l)2003l) . The resulting value 
is consistent with the observed hard X-ray luminosity of 
NGC 4435. Moreover, we notice that the X-ray lumi- 
nosity is much less dependent on the age of the stellar 
population (the turn-off mass of the stars) than the FIR 
luminosity. This implies that, in post starburst galaxies, 
the X-ray - FIR luminosity ratio tends to be larger than 
that of pure starbursts which are dominated by massive 
stars. 

Another indication of the strength of the AGN comes 
from the comparison of the hard X-ray luminosity and 
Ha emission. Using the relati on between Lx-Ln a , de- 
rived for AGN (jHo et alJl200l|) . the expected Ha lumi- 
nosity would be ^hS N ~ 6.6 x 10 31 W. This value is about 
15% of the observed Ln a corrected for extinction via the 
Balmer decrement. This contribution to the ionizing flux 
is higher than the limit derived from mid infrared nebu- 
lar lines (Sect. I3.2.1JI . However it must be considered as 
a generous upper limit because on one hand all the X-ray 
luminosity has been attributed to the AGN and, on the 
other, it is known that the Balmer decr ement provides 
only a lower limit to the attenuation fe.g. lPanuzzo et al.l 
2003). In fact if we correct the observed Ln a with the 
attenuation derived from our best fit model we obtain 
an upper limit to the AGN contribution ^5%, consistent 
with what is found from the [Ne in]/[Ne n] ratio. 



5.7. NGC 4-^35 and the phenomenology of the 
galaxy-galaxy interaction 

NGC 4435 is inte racting with NGC 44 38; in partic- 
ular, simulations by iVollmer et al.1 l)2005l) suggest that 
NGC 4435 passed through the disk of NGC 4438 about 
100 Myr ago at a radial distance of 5-10 kpc. With such 
an impact parameter the interaction of the ISM of the 
two galaxies is unavoidable. 

The age of the starburst t^ derived from our fit is 
larger than the (dynamically estimated) time since the 
ISM-ISM collision (around 100 Myr), but it is consistent 
with the epoch of the onset of interaction. This suggests 
that the star formation was triggered by the int eraction 
(as o ften seen in galaxy-galaxy simulations, e.g. iBarnesI 
I2004D possibly producing a relevant activity also before 
the ISM-ISM collision. We cannot, however, draw strong 
conclusions on the relative importance of the tidal inter- 
action and ISM-ISM collision, because it is difficult to 
derive the detailed shape of the star formation history 
during the burst. 

Another phenomenon that may be triggered by 
galaxy-galaxy inter action is AGN activity (see e.g. 
Hopk ins et al.ll200q) . However the various lines of evi- 
dence shown in this work indicate a lack of current AGN 
activity. 

The age of the circumnuclear starburst is several times 
larger than the typical duration of the AGN phase and 
we cannot rule out the possibility that an AGN was once 
active in NGC 4435. However the surrounding gas is 
predominantly in a cold phase w ith only a fraction in a 
hot component (~ 3 x 10 7 M , iMachacek et all 120041) 
and it is supported by circular motion " ifCoccato et alJ 
2006), indicative of a relaxed ISM. We conclude that the 
feedback from the possible AGN was not strong enough 
to halt the star formation process which is now fading as 
a result of the gas consumption. 

Alternatively, the AGN activity has not been trig- 
gered yet (and perh aps will not be triggered at all). 
iCoccato et al.1 l)2006l) provide a 3tr upper mass limit of 
7.5xl0 6 M Q for the possible central SMBH, an order 
of magnitude lower than the value obtained from the 
Mbh-o- and AjfBH-^buigc relations. In this case NGC 4435 
could be a prototype of the so called "laggard" SMBH 
which are believed to undergo slow gas fueling of the 
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galactic center and limited growth of the SMBH mass 
l|Vittormi et al1l200l . 

The above picture sketched for the evolution of 
NGC 4435 would suggest that the star formation activ- 
ity could inhibit the AGN feeding both by consuming the 
gas and by SNe feedback. Both cases point towards a rel- 
atively low efficiency of the AGN feedback with respect 
to that of the SNe in i ntermediate /small m ass galaxies 
like NGC 4435 (see e.g. IGranato et alJl200^ . 

6. CONCLUSIONS 

We have performed a thorough spectrophotometric 
study of the early-type galaxy NGC 4435 using new 
Spitzer IRS spectra combined with IRAC and MIPS 
archival data and existing broad band measurements 
from X-ray to radio wavelengths. 

The IRS spectrum is dominated by PAH features, neb- 
ular emission lines and molecular hydrogen rotational 
lines that originate in the circumnuclear disk that char- 
acterizes this galaxy. 

The PAH features in the spectrum are very similar to 
those observed in other star forming galaxies. In particu- 
lar, we repor t the detection of the PAH feature observed 
at 17.1 /xm (|Smith et al.ll2004l) which, due to the lower 
contamination by the H2 S(l) rotational line, we show is 
actually a blend of two distinct components. 

The panchromatic SED of the central 5" region was 
analysed with the spectrophotometric code GRASIL. 

This SED is well reproduced from the UV to radio by 
a model of a fading starburst superimposed on a simple 
stellar population of 8 Gyr and solar metallicity. 

The analysis of MIR nebular lines, combined with that 
of the SED, allow a direct estimate of the gas metallicity, 
which also results as solar. 

Though a low level nucle ar activity has been su spected 
in N GC 4435 from optical l)Ho et al.lll997ll2002|) and X- 
ray l|Machace k et al .1120041) observations, we collect the 
following evidence for the lack of AGN excitation: 

(i) we have not detected high excitation nebular emis- 
sion lines in the MIR spectrum 

(ii) the [Neni]15.5/[Neii]12.8 ratio constrains the con- 
tribution of a possible AGN to the ionizing flux to be less 
than 2%. 

(Hi) the upper limit on the temperature derived from 
H2 S(l) and S(2) rotational lines is lower than expected 
for AGN excitation. 



(iv) the X-ray emission is within the range expected 
from X-ray binaries from the starburst 

The age of the starburst, 186 ±37 Myr, corresponds to 
the epoch of the last encounter with NGC 4438 derived 
from dynamical simulations, suggesting that this inter- 
action has triggered the observed star formation event 
ijCombes et al.lll988|) . However we cannot draw strong 
conclusions on the relative importance of the tidal inter- 
action and ISM-ISM collision, because of the difficulty in 
deriving the detailed shape of the star formation history 
during the burst. 

The combination of Spitzer observations with the anal- 
ysis of the complete SED allows an accurate determina- 
tion of the mass and epoch of the rejuvenation episode. 
The mass of the starburst (~ 1.40 x 10 8 M Q ) corresponds 
to about 1.4% of the stellar mass sampled by the central 
5 arcsec aperture. The young component makes the op- 
tical spectrum of NGC 4435 closely similar to a typical 
intera cting early- type galaxy with inverted C a n[H+K] 
lines l)Lon£fhetti et al]ll999t iSarzi et al.ll2005[) that will 
later turn into a typical cluster E+A galaxy. This study 
supports the notion that early-type galaxies with rela- 
tively strong hydrogen absorption features are due to re- 
cent small rejuvenation episodes, r ather than being th e 
result of delayed galaxy formation ijBressan et alj fl9961. 

Finally the total SED of NGC 4435 is fitted by simply 
adding to the starburst a larger fraction of the old stel- 
lar population. The starburst provides only 0.3% of the 
total stellar mass of the galaxy but 8% of the bolometric 
luminosity. 
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